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Abstract 

We report the observation of a narrow charmonium-like state produced in the exclusive decay 
process K^7r~^7r~ J/tjj. This state, which decays into tt'^tt" J/V', has a mass of 3872.0 ± 

0.6(stat) lb 0.5(syst) MeV, a value that is very near the Md + Mn* mass threshold. The results are 
based on an analysis of 152M B-B events collected at the T(4S') resonance in the Belle detector 
at the KEKB collider. The signal has a statistical significance that is in excess of 10a. 

PACS numbers: 14.40.Gx,12.39.Mk,13.20.He 
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A major experimental issue for the cc charmonium particle system is the existence of as 
yet unestablished charmonium states that are expected to be below threshold for decays to 
open charm and, thus, narrow. These include the n = 1 singlet P state, the J^^ = l^Pci, 
and possibly the n = 1 singlet and triplet spin-2 D states, i.e. the J^^ = 2"+ 1^D^2 and 
J^*^ = 2 l^Dc2, all of which are narrow if their masses are below the DD* threshold. 
The observation of these states and the determination of their masses would provide useful 
information about the spin dependence of the charmonium potential. 

In addition to charmonium states, some authors have predicted the existence of 
"molecular charmonium" states [l] and ccg "hybrid charmonium" states 2]. If such states 
exist with masses below the relevant open charm threshold, they are expected to be narrow 
and to have large branching fractions to low-lying cc charmonium states. 

The large B meson samples produced at i?-factories provide excellent opportunities to 
search for new charmonium states. The Belle group recently reported the first observation 
of the ric{2S) via its KsK~tt~^ decay channel in exclusive B — ^ KKsK~tt~^ decays based 
on an analysis of 44. 8M B-B events 0. Strategies for finding the remaining missing states 
have been presented by Eichten, Lane and Quigg 0]; they note that a narrow ^Dc2 should 
have substantial decay branching fractions for 7Xci and tt+tt" J/'0 final states. In this paper, 
we report on an experimental study of the h^ti' J/ip and 7Xci niass spectra from exclusive 
B^ K^TT^TT" J/ip and K'^'yxci decays ^ using a 152M B-B event sample. The data 
were collected in the Belle detector at the KEKB energy-asymmetric e^e~ collider, which 
operates at a center-of-mass (CM) energy of a/s = 10.58 GeV, corresponding to the mass of 
the T(4S') resonance. KEKB is described in detail in ref. @|. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector, a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOP), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect Kl mesons and to identify 
muons (KLM). The detector is described in detail elsewhere 0. 

The identity of each charged track is determined by a sequence of likelihood ratios that 
determine the hypothesis that best matches the available information. Tracks are identified 
as pions or kaons based on the specific ionization in the CDC as well as the TOP and ACC 
responses. This classification is superseded if the track is identified as a lepton: electrons 
are identified by the presence of a matching ECL cluster with energy and transverse profile 
consistent with an electromagnetic shower; muons are identified by their range and transverse 
scattering in the KLM. 

Por the B — » Ktt^tt~J/iP study we use events that have a pair of well identified oppositely 
charged electrons or muons with an invariant mass in the range 3.077 < M£+£- < 3.117 GeV, 
a loosely identified charged kaon and a pair of oppositely charged pions. In order to reject 
background from 7 conversion products and curling tracks, we require the vr+vr" invariant 
mass to be greater than 0.4 GeV. To reduce the level of e"''e~ qq {q = u, d, s or c-quark) 
continuum events in the sample, we also require R2 < 0.4, where R2 is the normalized 
Pox-Wolfram moment |^, and | cos9b\ < 0.8, where 9b is the polar angle of the 5-meson 
direction in the CM frame. 

Candidate B^ — > K^ti^'k^ J/ip mesons are reconstructed using the energy difference 

AE = E'j^^ - -EbiLi and the beam-energy constrained mass Mbc = \/ {E^^^)"^ - {p'^Y, 
where -Ebeam is ^^e beam energy in the CM system, and E'^ and are the CM energy 
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and momentum of the B candidate. The signal region is defined as 5.271 GeV < Mbc < 
5.289 GeV and \AE\ < 0.030 GeV. 



^ 12000 




0.40 0.80 1.20 0.40 0.80 1.20 



M(jtvrr) - M(rr) (GeV) M(jivrr) - M(rr) (GeV) 

FIG. 1: Distribution of M(7r+7r~£+^^) — M{£~^i~) for selected events in the AE-M^^ signal region 
for (a) Belle data and (b) generic B-B MC events . 

Figure HJ^a) shows the distribution of AM = M{7r~^7i~i~^£~) — M{i~^£~) for events in 
the AE-M\yc signal region. Here a large peak corresponding to ip' —>■ 'K^Ti^J/ip is evident 
at 0.589 GeV. In addition there is a significant spike in the distribution at 0.775 GeV. 
Figure E^b) shows the same distribution for a large sample of generic B-B Monte Carlo 
(MC) events. Except for the prominent ip' peak, the distribution is smooth and featureless. 
In the rest of this paper we use M^Tt^-K' J determined from AM + Mjj^, where Mj/^ is 
the PDG [9] value for the J /if) mass. The spike at AM = 0.775 GeV corresponds to a mass 
near 3872 MeV. 

We make separate fits to the data in the ip' (3580 MeV < M.^+^-j/^ < 3780 MeV) and 
the M = 3872 MeV (3770 MeV < M^+^~j/^ < 3970 MeV) regions using a simultaneous 
unbinned maximum likelihood fit to the Mbc, AE, and Ai^+^-j/^ distributions For 
the fits, the probability density functions (PDFs) for the Mbc and M.^+^-j/.0 signals are 
single Gaussians; the AE signal PDF is a double Gaussian comprised of a narrow "core" 
and a broad "tail" [l7]. The background PDFs for AE and M^+i^-jj^ are linear functions, 
the Mbc background PDF is the ARGUS threshold function For the ip' region fit, the 
peak positions and widths of the three signal PDFs, the AE core fraction, as well as the 
parameters of the background PDFs are left as free parameters. The values of the resolution 
parameters that are returned by the fit, listed in Table HJ are consistent with MC-based 
expectations. For the fit to the M = 3872 MeV region, the Mbc peak and width, as well as 
the AE peak, widths and core fraction are fixed at the values determined from the ip' fit- 

The results of the fits are presented in Table ITTl Figures |2f a), (b) and (c) show the 
Mbc, M.^+^-j/^, and AE signal-band projections for the M = 3872 MeV signal region, 
respectively. The superimposed curves indicate the results of the fit. There are clear peaks 
with consistent yields in all three quantities. The signal yield of 35.7 ± 6.8 events has a 

statistical significance of 10.3(j, determined from ^—2 ln(£o/'^max), where £max and Cq are 
the likelihood values for the best-fit and for zero-signal-yield, respectively. In the following 
we refer to this as the X(3872). 
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TABLE I: Resolution values from the fits to the ip' signal region. The errors are statistical only. 



Quantity 


Fitted value 




2.6 ±0.1 MeV 


(TAB (core) 


11.6 ±0.4 MeV 


(TAB (tail) 


130 ± 130 MeV 


Core fraction 


0.965 ±0.015 




5.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 5.3 
Mbc (GeV) 




3.82 3.84 3.86 3.88 3.9 3.92 

M(J/v|; kk) (GeV) 




0.05 0.1 0.15 0.2 
AE (GeV) 



FIG. 2: Signal-band projections of (a) Mbc, (b) M^+^-j/^ and (c) AE for the X(3872) 
7r~^ iT~ J/ ip signal region with the results of the unbinned fit superimposed. 



We determine the mass of the signal peak relative to the well measured ip' mass: 



^meas _ ^meas ^ ^PDG _ 3372.0 ± 0.6 ± 0.5 MeV. 



Here the first error is statistical and the second systematic. Since we use the precisely known 
value of the ip' mass as a reference, the systematic error is small. The M^/ measurement, 
which is referenced to the J/ip mass that is 589 MeV away, is — 0.5±0.2 MeV from its world- 
average value Variation of the mass scale from M^/ to Mx requires an extrapolation 
of only 186 MeV and, thus, can safely be expected to be less than this amount. We assign 
0.5 MeV as the systematic error on the mass. 

The measured width of the X(3872) peak is (T = 2.5 ± 0.5 MeV, which is consistent 
with the MC-determined resolution and the value obtained from the fit to the ip' signal. 
To determine an upper limit on the total width, we repeated the fits using a resolution- 



TABLE II: Results of the fits to the ip' and M = 3872 MeV regions. The errors are statistical only. 



Quantity 


Ip' region 


M = 3872 MeV region 


Signal events 


489 ± 23 


35.7 ±6.8 


^X-j/^ peak 


3685.5 ± 0.2 MeV 


3871.5 ± 0.6 MeV 


^MlT+TT-J/lp 


3.3 ± 0.2 MeV 


2.5 ± 0.5 MeV 
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broadened Breit-Wigner (BW) function to represent the signal. This fit gives a BW width 
parameter that is consistent with zero: F = 1.4 ± 0.7 MeV. From this we infer a 90% 
confidence level (CL) upper limit of F < 2.3 MeV. 

The open histogram in Fig. |21^a) shows the vr'^vr" invariant mass distribution for events in 
a ±5 MeV window around the X(3872) peak; the shaded histogram shows the corresponding 
distribution for events in the non-signal Ai?-Mbc region, normalized to the signal area. The 
TT+TT" invariant masses tend to cluster near the kinematic boundary, which is around the p 
mass; the entries below the p are consistent with background. For comparison, we show the 
TT+TT^ mass distribution for the tp' events in Fig. Efb), where the horizontal scale is shifted 
and expanded to account for the different kinematically allowed region. This distribution 
also peaks near the upper kinematic limit, which in this case is near 590 MeV. 



o) X(3S72) 




0.41 0.51 
M(7t*7i ) (GeV) 



FIG. 3: Af(7r^7r^) distribution for events in the (a) M [tt'^ it~ J / ip) = 3872 MeV signal region, and 
(b) the ip' region. The shaded histograms are sideband data normalized to the signal-box area. 
Note the different horizontal scales. 

We determine a ratio of product branching fractions for i^'+X(3872), X(3872) — »• 

Tx^-K^ J/ip and —>■ K^ip', ip' Ti^'K^J/tp to be 

B{B+ ^ K+X (3872)) X B{X (3872) -^71+ 71- J/ ij) ^ , n ^ 

^ — — V^, /! f — - — — = 0.063 ± 0.012(stat) ± 0.007(syst). 

Here the systematic error is mainly due to the uncertainties in the efficiency for the 
X(3872) 7c^7[-J/ip channel, which is estimated with MC simulations that use differ- 
ent models for the decay . 



The decay of the ^Dc2 charmonium state to 7Xci is an allowed El transition with a 
partial width that is expected to be substantially larger than that for the 7i^7i~J/ip final 
state; e.g. the authors of Ref. 0] predict F(^Dc2 IXa) > 5 x F(^Dc2 tt+tt" J/-?/^). 
Thus, a measurement of the width for this decay channel can provide important information 
about the nature of the observed state. We searched for an X(3872) signal in the 7Xci decay 
channel, concentrating on the Xci iJ 1'^ final state. 

We select events with the same J/ip ^ t^t~ and charged kaon requirements plus two 
photons, each with energy more than 40 MeV. We reject photons that form a tt*^ when 
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combined with any other photon in the event. We require one of the j'^ combinations to 
satisfy 398 MeV < {M^i^i- - M^+^-) < 423 MeV (corresponding to -15 MeV < (M^j/^ 



M^^^) < 10 MeV). In the following we use M 



7Xcl 



■M. 



+ MPp^, where M™g 



Xcl 



is the PDG Xci mass value 

The B K'yXci, Xci iJ/i' decay processes have a large combinatoric background 
from B — > Kxci decays plus an uncorrelated 7 from the accompanying B meson. This 
background, which tends to peak in the Mbc signal region, produces a peaking in AE at 
positive values that is well separated from zero, and is removed by the AE < 30 MeV 
requirement. Because of the complicated AE background shape and its correlation with 
Mbc, we do not include AE in the unbinned likelihood fit. Instead we do an unbinned fit 
to the M^^^-^^ and Mbc distributions with the same signal and background PDFs for Mbc 
and M^^^-^ that are used for the n^Tc^J/ip fits. We fix the values of the Gaussian widths at 
their MC values, and the ip' and X(3872) masses at the values found from the fits to the 
71^71^ J/ip channels. The signal yields and background parameters are allowed to float. 

The signal-band projections of Mbc and M^^^-^ for the ip' region are shown in Figs. IH^a) 
and (b), respectively, together with curves that represent the results of the fit. The fitted 
signal yield is 34.1 ± 6.9 ±4.1 events, where the first error is statistical and the second is 
a systematic error determined by varying the Mbc and M^-^^-^ resolutions over their allowed 
range of values. The number of observed events is consistent with the expected yield of 
26 ± 4 events based on the known B Kip' and ip' IXa branching fractions JP] and the 
MC-determined acceptance. 




i.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 5.3 

Mbc (GeV) 




3.58 3.6 3.62 3.64 3.66 3.68 3.7 3.72 3.74 3.76 3.78 

M(^,) (GeV) 




5.2 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 5.3 

Mbc (GeV) 




3.78 3.8 3.82 3.84 3.86 3.88 3.9 3.92 3.94 3.9- 

M(l7cl) (GeV) 



FIG. 4: Signal-band projections of (a) Afbc and (b) M^^^^ for the tp' region with the results of 
the unbinned fit superimposed, (c) and (d) are the corresponding results for the M = 3872 MeV 
region. 



The results of the application of the same procedure to the M = 3872 MeV region are 
shown in Figs. |^c) and (d). Here, no signal is evident; the fitted signal yield is 3.7 ± 3.7 ± 
2.2 events, where the first error is statistical and the second is a systematic error determined 
by varying the input parameters to the fit over their allowed range. From these results, we 
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determine a 90% CL upper limit on the ratio of partial widths of 



r(X(3872) ^ TT + TT- J/^/-) 



where the effects of systematic errors have been included. This limit on the 7Xci decay 
width contradicts expectations for the ^T)c2 charmonium state. 

The mass of the observed state is higher than potential model expectations for the center- 
of-gravity (cog) of the l^Dcj states: the Cornell 1^ and the Buchmiiller-Tye [l^ potentials 
both give Mm^dD) = 3810 MeV, which is 60 MeV below our measurement. A model by 
Godfrey and Isgur |l2| gives a ^Dc2 mass of 3840 MeV but also predicts a ^Dd of 3820 MeV, 
which is higher than its observed value of 3770 MeV. Identification of the observed state 
with the ^Dc2 state would imply a ^Dc2 — ^ Dd splitting of ~ 100 MeV. On the other hand, 
complications due to effects of coupling to real DD and virtual DD* states may reduce the 
validity of potential model calculations for these states [iS^ . 

The 7r"^7r~ invariant masses for the M = 3872 MeV signal region concentrate near the 
upper kinematic boundary as is also the case for ip' -n^n^ J /ip (see Figs. Efa) and (b)). 
For the M = 3872 MeV signal region, however, the boundary corresponds to the mass of 
the p. The decay of a cc charmonium state to p J/ip is an isospin- violating process, and 
these are strongly suppressed in ip' ^ J/ip transitions. With more data, we will be able to 
determine whether or not the tt+tt^ system is coming from p meson decay jl^ and thereby 
establish the C parity of the X(3872). Information about other possible decay channels, 
such as DD^DDtt^ and DD'y, would be useful for determining other quantum numbers of 
this state jl9|. 

In summary, we have observed a strong signal for a state that decays to -n-^n^ J/ip with 

M = 3872.0 ± 0.6 (stat) ± 0.5 (syst) MeV 
r < 2.3 MeV (90% CL). 

This mass value and the absence of a strong signal in the 7Xci decay channel are in some 
disagreement with potential model expectations for the ^Dc2 charmonium state. The mass 
is within errors of the D^D*^ mass threshold (3871.1 ± 1.0 MeV l^j), which is suggestive of 
a loosely bound DD* multiquark "molecular state," as proposed by some authors [L ,2(3|. 
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